Time-resolved imaging of magnetoelectric domain switching in multiferroic M11WO4 
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We show that the electric-field-induced reversal of the magnetic order parameter in multiferroic MnW04 
occurs on the time scale of milliseconds. Throughout the switching process the magnetic order and the magnet- 
ically induced electric polarization remain rigidly coupled. The temporal progression of the domain structure 
was imaged with nanosecond resolution by an electrical-pump-optical-probe technique using optical second 
harmonic generation. These nonequilibrium domain states significantly differ from the quasi-static domain 
reversal. A qualitative model gives an estimate of why the magnetoelectric order-parameter reversal in the 
magnetically induced ferroelectrics is not inherently ultrafast. 

PACS numbers: 75.85.+t, 75.78.Fg, 77.80.Fm, 42.65.Ky 



I. INTRODUCTION: TIME-RESOLVED MULTIFERROICS 

Materials with cross-coupled magnetic and electric proper- 
ties, called magnetoelectrics, are intensely discussed because 
of their potential for controlling dielectric or magnetic prop- 
erties with the crosswise magnetic or electric field ? Magne- 
toelectric effects are intrinsically strong in compounds where 
magnetic long-range order breaks the inversion symmetry and 
induces a spontaneous polarization. This mechanism results 
in a rigid coupling of the improper ferroelectric polarization 
to the proper magnetic order parameter(s). 3 5 

Despite intensive research on these so-called "multifer- 
roics" one central aspect of the interplay of magnetic and elec- 
tric order parameters has scarcely been discussed: How does 
the voltage-induced reversal of the magnetization progress 
with time? Knowing the dynamics of voltage-induced magne- 
tization reversal would highlight the competition of spins and 
charges. This would significantly advance our general under- 
standing of the parameters defining the magnetoelectric cou- 
pling in one of the most prominent classes of multiferroic s. 
But also with respect to possible applications like memory 
devices and sensors the switching dynamics is a key issue. 
For example, the feasibility of any candidate compound for 
a novel type of memory element depends crucially on its bit- 
reversal time. 

Thus far, investigations of the dynamical properties of mul- 
tiferroics are devoted to the AC magnetoelectric response 
of magnetostrictive-piezoelectric composites at microwave 
frequencies. 6 In two multiferroics the ultrafast response to 
laser pulse excitation has been investigated.™ Yet, in none 
of these experiments time-resolved magnetoelectric switching 
processes were observed or discussed. In spite of the signifi- 
cance of this topic investigations are virtually non-existent. 

In this Report we show that the electric-field-induced re- 
versal of the magnetic order in multiferroic MnWC>4 extends 
over a surprisingly long time span in the order of milliseconds. 
The evolution of the multiferroic domain structure throughout 
the switching process is spatially resolved by optical second 
harmonic generation and compared to the quasi- static reversal 
of domains. Fundamental differences between the observed 
magnetoelectric switching and models for switching of do- 



mains in non-multiferroics are highlighted. 

II. SAMPLES AND METHODS 
A. Multiferroic MnW0 4 

The compound investigated here, MnWC>4, is prototypi- 
cal for the group of magnetically induced ferroelectric sPESl 
There is a multitude of similar compounds promoting non- 
centro symmetric, preferably frustrated and/or incommensu- 
rate spin structures. They can induce spontaneous polariza- 
tion even up to ambient temperature!^ This background em- 
phasizes the general merit of our results. 

Three magnetically ordered phases discussed in detail in 
Refs. 12l and 1131 constitute the phase diagram of MnWCU. In 
the AF3 phase below 7]y = 13.5 K the magnetic moments 
of Mn 2+ align collinearly along the easy axis while their 
magnitudes are sinusoidally modulated yielding incommen- 
surate order in the xz plane (axes: see Ref. [331). An addi- 
tional transverse spin component orders at T2 = 12.7 K. In 
this AF3 —> AF2 transition the spin-density wave turns into 
an incommensurate elliptical spin spiral that breaks the in- 
version symmetry and induces a spontaneous polarization Py P 
along the y axis. Another transition at T\ =7.6 K leads to 
the AF1 phase with a collinear centrosymmetric alignment of 
spins along the easy axis in the xz plane. The multiferroic AF2 
phase is characterized by two magnetic domain states with op- 
posite helicity of the associated spin spiral. The helicity o± is 
coupled to the direction of the spontaneous polarization ±Py V 
so that both can be set or reversed by an electric field ±E y ¥^ 

B. Optical second harmonic generation 

A convenient way for imaging the multiferroic domain 
structure is optical second harmonic generation (SHG). Laser 
light is incident onto a crystal where it induces a frequency 
doubled wave (Fig.[TJe)). The polarization of the SHG wave 
is determined by the symmetry of the crystal and, hence, by 
any form of long-range order affecting the symmetry. Thus, 
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SHG suggests itself for probing (anti-) ferroic order of any 
type, especially in compounds with multiple order parame- 
ters. In particular, the optical approach allows one to in- 
vestigate the ordered state with high temporal and spatial 
resolution 6 16 and resolve the magnetoelectric switching of 
the domains. In previous SHG experiments, the as-grown 
multiferroic domain structure in MnWCU was investigated 17 
and the electric-field-induced conversion into a single-domain 
state was observed in pyroelectric, neutron diffraction, and 
SHG experiments ! 1 5 1 1 8 1 1 9 1 However, in none of these experi- 
ments the actual poling process with its temporally evolving 
domain structure was resolved. 



C. Experimental details 

The electric field was applied along the y axis of the sam- 
ples. As shown in the insets of Fig.[2| this is achieved for the 
(010) samples by glass plates coated with 100 nm of transpar- 
ent indium-tin-oxide, whereas for the (100) samples polished 
steel electrodes are used. The quasi- static evolution of the do- 
main structure in Fig. [T] was investigated by applying the elec- 
tric field to the electrodes with a tunable high-voltage supply. 
The dynamic reversal of the domains in Figs. [2] to |4] was re- 
solved with the electrical-pump-optical-probe technique de- 
picted in Fig. |2ja). Repetitive rectangular voltage pulses are 
applied to the samples. At time zero the voltage pulse re- 
verses the electric field applied to the sample from —E y to 
+E y as indicated by the bar on top of Figs. [2|a), |2jb), and 
[2jc). The slope of the voltage pulse is 50 ns. Different fields 
up to E y = ±900 kV/m were used. The state of the sample is 
probed by SHG at the time At while the electric field pulse is 
still applied. After 50 ms the electric field is reversed to —E y 
again. However, this exceeds the time span displayed in our 
figures so that we only see the initial reversal of the electric 
field at t = 0. SHG data are accumulated across about 1000 
laser pulses and, thus, voltage cycles per data point. 

The samples are placed in an optical cryostat. The tem- 
poral evolution of the magnetic system is tracked by chang- 
ing the delay At between the voltage reversal at t = and 
the 5-ns laser pulses (energy: 1—5 mJ) as source of the 
SHG. The transmission setup for magnetically induced SHG 
and a comprehensive characterization of SHG on multifer- 
roic MnWCU are detailed in Refs. and |20| Summarizing 
these, the SHG experiments were done on MnWO^OlO) sam- 
ples in which the X^ z > contribution at 2.21 eV couples to 
the incommensurate magnetic order but not to the sponta- 
neous polarization. 14 20 In addition, the magnetically induced 
polarization was probed via the contribution at 2.72 eV 
in MnWO4(100) samples. 20 All samples were grown by top- 
seeded melt growth and lapped and etch-polished to a thick- 
ness of 0.4 — 1.0 mm with a silica slurry. A summary of the 
samples, SHG contributions, and types of order probed in our 
experiments is given in Table [I] 



III. RESULTS AND DISCUSSION 
A. Quasi-static switching of domains 

Figure [T] shows the quasi-static switching of magnetic he- 
licity domains in MnWO4(010) by an electric field. Starting 
from a zero-field-cooled state (E y = 0) or a single-domain 
state (E y = —750 kV/m), the electric field in Figs. [TJa) to 
[T|c) and Figs.[TJf) toJTJh), respectively, is increased to E y = 
+750 kV/m with a time > 1 s between subsequent data points. 
The corresponding virgin curve in Fig. [TJd) is followed by a 
hysteresis loop that is well centered in contrast to the loops 
in Ref. [21] The coercive field is 120 kV/m and the largest 
field applied exceeds the saturation field of ±600 kV/m by 
25%. The domain structure obtained after zero-field cooling 
was already discussed in Ref. [TTJ A magnetic bubble topol- 
ogy emerges. On the surface of the (010) sample the domain 
boundaries yield an affinity to be oriented along the z axis and 
the magnetic easy axis. Domain states of opposite helicity dif- 
fer in their brightness because of SHG interference effects .™ 
In the present work we see that in the increasing electric field 
the bubbles for either of the domain states shrink with a ten- 
dency to minimize the surface. This leads to the spherical 
domain shape visible on the (010) face in Fig.[TJc). Finally, a 
single-domain state is acquired. 

In stark contrast to the bubble topology the re-emergence 
of the opposite domain state during gradual reversal of the 
poling field leads to the formation of domains that appear as 
stripes along the z axis on the surface of the (010) sample, 
see Figs.[TJf) to[TJh). Nucleation of domains occurs primarily 
at the sample edges but was also observed within the sample. 
Note that the direction of the stripes is not that of the sponta- 
neous polarization or of the electric poling field. 

As an interpretation, the nucleation of domains in Figs.[TJa) 
to [TJc) is determined by the emergence of a spin spiral upon 
crossing T2. This is an entirely magnetic phenomenon be- 
cause the primary order parameter is magnetic and electric 
fields are not involved. In contrast, the domain structure in 
Figs.[TJf) to[TJh) results from the modification of a spin spiral 
that is already present. This modification is based on magne- 
toelectric interactions because it is driven by the electric field. 
Note that the weakest magnetic exchange is found along the x 
direction!^ This explains the formation of the domain stripes 
along z on the (010) face. 

B. Dynamic switching of domains 

Figure [2] shows the dynamic switching of magnetic helicity 
domains in MnWO4(010) by an electric field. As mentioned 
in Section HTCl the fastest time in which the electric field can 
be reversed is 50 ns. This is determined by the slope of the 
repetitively applied voltage pulses. Figure |2jb) reveals that 
the response of the magnetic order parameter to the electric 
field pulses of E y = 750 kV/m is strikingly slow: The rever- 
sal of the magnetic domain state takes as long as 20 ms. The 
domain structure throughout this dynamic reversal is shown 
in Figs. [2fd) to|2jf). It differs substantially from the one in 
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FIG. 1 : (Color online) Evolution of the magnetic domain structure 
in multiferroic MnWO^OlO) during quasi-static electric-field pol- 
ing, (a-c) Poling from a multidomain state after zero-field cool- 
ing, (d) SHG hysteresis loop with open symbols denoting the vir- 
gin curve, (e) Sample environment with applied voltage (HV) and 
incident fundamental and emitted SHG light, (f-h) Poling from a 
single-domain state at —750 kV/m during gradual increase of the 
field towards +750 kV/m. Bright and dark areas in the SHG images 
correspond to domains with opposite magnetic helicity. Temperature 
isllK. 
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FIG. 2: (Color online) Evolution of the domain structure during dy- 
namic electric-field poling with a voltage pulse (E = 750 kV/m, slope 
50 ns, width 50 ms). (a) Sketch of the repetitive electric-field poling 
with subsequent SHG probing, (b, c) SHG intensity as a function of 
the delay to the voltage pulse at t = for (b) MnWO 4 (010) 11.8 K, 
750 kV/m and (c) MnWO 4 (100) 12.0 K, 330 kV/m. Insets show the 
position of the electrodes with respect to the laser beam. The dashed 
lines show the SHG intensity for the two single domain states. Note 
that the abrupt jump in (c) at t = is caused by EFISH and is thus 
not related to the long range order, (d-f) Magnetic domain structure 
at various delays. Arrows: see text. Bars on top of (a-c) indicate the 
applied electric field. 



zero-field-cooled samples and also from the one obtained dur- 
ing quasi- static field reversal (Fig. [TJ. The dynamic reversal 
leads to the nucleation of rhomboid- shaped domains whose 
sides run approximately parallel to the z axis and the mag- 
netic easy axis. In this nonequilibrium state, a tendency to 
minimize the domain surface is already present in the form of 
rounded corners (arrow in Fig. [2jf» and domain fusion (ar- 
row in Fig.[2je)) but much less pronounced than in Fig.[TJc). 
Although the images are accumulated from about 1000 laser 
pulses (i.e. 1000 poling cycles) a well-defined domain struc- 
ture is obtained. Therefore, the domain reversal always pro- 
gresses through the same domain pattern. This pronounced 
memory effect is a strong indication for pinning effects. 

The same type of dynamic switching experiment was also 
carried out on a MnWO4(100) sample with the result shown 
in Fig. [2j c )« m this geometry the electric field is limited to 
a maximum value of ±330 kV/m because of the larger dis- 
tance between the electrodes and the limited voltage of the 



HV source. At t < the SHG intensity remains constant for 
an electric field of —330 kV/m. The electric-field reversal to 
330 kV/m at t = leads to a jump of the SHG signal within 
the rise time of the voltage pulse. This is followed by a grad- 
ual decrease of the SHG intensity lasting until t = 10 ms. At 
t > 10 ms, the SHG intensity remains constant at a level that 
is different from that at t < 0. Upon the field reversal from 
+330 to —330 kV/m (not shown) the SHG intensity returns to 
its original value. 

The jump at t = is caused by SHG contributions that 
are induced by the externally applied field but unrelated to 
the long-range order. 19 This effect is well-known and docu- 
mented as "EFISH" (electric-field-induced second harmonic) 
in the literature. 23 The EFISH signal follows the applied volt- 
age without delay so that its rise time is 50 ns in Fig. |5J C )- 
EFISH contributions are always polarized parallel to the ap- 
plied electric field which explains their absence in Fig. |2|b) 
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where the light is propagating along the field direction. Only 
the decrease of the SHG signal from the upper to the lower 
dashed line in Fig.[2jc) is associated to the actual domain re- 
versal. The reversal changes the domain-related SHG con- 
tribution. This changes the interference with the (constant) 
EFISH contribution and, hence, the net SHG yield until it sta- 
bilizes after the domain reversal is completed at t = 10 ms. 
The factor of two in switching time between Figs. |2jb) and 
[2jc) originates from differences in the sample temperature and 
the applied electric fields, as discussed in detail below. 

As mentioned, the (100) sample promotes SHG contribu- 
tions coupling to the spontaneous polarization P3 Any tran- 
sient decoupling of the polarization from the magnetization 
would therefore manifest in Fig. |2jc) as a timescale differ- 
ent from that of the magnetization in Fig.[2jb). In particular, 
a fast response was considered likely for any electronic con- 
tributions to the spontaneous polarization ! 24 ! 25 ! However, no 
such decoupling is observed. (Remember that the time scale 
introduced by the jump at t = is related to the applied elec- 
tric field and not to a switching of the multiferroic order.) We 
thus conclude that the magnetic and electric order parameter 
remain rigidly clamped throughout the switching process so 
that the polarization, just like the magnetization, reverses on 
a millisecond time scale. Since faster switching processes are 
not present it is not necessary to increase the time resolution of 
the experiment beyond the nanosecond range employed here. 

In order to investigate variations of the magnetoelectric 
domain-reversal time t dynamic switching experiments as in 
Fig. [2jb) we carried out at different temperatures and electric 
fields. The accessible temperature interval is limited to 0.6 K 
by the increase of the coercive field towards lower and the de- 
crease of the SHG intensity towards higher temperature. The 
electric field is determined by the field required to acquire a 
single-domain state and is limited to 900 kV/m by the break- 
through field in helium gas. Figure [3] shows the reversal time 
in dependence of these parameters, probed via the magnetic 
SHG signal on the (010) sample. First, Fig.[3ja) reveals that T 
increases with decreasing temperature. Although a tempera- 
ture interval of only 0.6 K is investigated T changes by a factor 
of 40, from 1 ms at 12.0 K to 40 ms at 11.4 K in an electric 
field of 750 kV/m. Figure [3jb) shows a series of experiments 
at 12.0 K with fields between 300 kV/m (the saturation field 
at 12.0 K) and 900 kV/m. Clearly, t decreases with increas- 
ing field. However, due to the limited field range we cannot 
make a conclusive statement on the validity of Merz's 
which describes the field dependence of the switching time in 
ferroelectrics. It predicts a relation % oc exp(—E/E a ) with E a 
as activation field. 

In Fig. [4] the evolution of two selected domains in a 
MnWO4(100) sample for a scenario with two field reversals at 
t = and at t = 20 ms is investigated. The area of the domains 
is shown in dependence of the delay after the first field rever- 
sal; the respective electric field is indicated by the bar on top 
of Fig.Qa). The timing of the voltage pulses was chosen such 
that the second field reversal occurs before the first domain 
reversal has been completed. After the initial voltage rever- 
sal at t = the areas of domains 1 and 2 grow steadily, albeit 
at different rates. After the second voltage reversal, domain 
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FIG. 3: Dependence of the magnetoelectric domain-reversal time on 
(a) temperature and (b) electric field. 



1 shrinks again, with a small discontinuity right at the volt- 
age reversal. In contrast, domain 2 disappears abruptly, and 
unrelated to the voltage reversal. Hence, although the mag- 
netoelectric domain switching progresses smoothly according 
to Fig. [2] Fig. [4] reveals that local deviations from the aver- 
aged response are possible. However, these local deviations 
are not predictable or controllable and they change with each 
annealing cycle applied to the sample. Note that the collapse 
of domain 2 shows that the intrinsic speed of domain switch- 
ing is orders of magnitude faster than the average response 
of the domains to the electric field. Collapsing domains were 
observed on the (010) as well as on the (100) face of the sam- 
ple with collapse times faster than 25 ns, the smallest delay 
between two images investigated in this experiment. 
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FIG. 4: Evolution of the domain structure in MnWOzKlOO) during 
repeated electric-field poling with E = ±330 kV/m at 11.5 K. (a, b) 
Area of domains 1 and 2 as a function of the delay to the first voltage 
pulse. The bar on top indicates the applied electric field, (c-e) SHG 
images of the domains at the given delays. 
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IV. DISCUSSION OF THE MAGNETOELECTRIC 
SWITCHING 

Summarizing Figs. [T] to [4] we conclude that the reversal of 
magnetic domains by an electric field is governed by domain 
wall movements on a timescale in the order of milliseconds 
and thus a surprisingly slow process. Pronounced memory ef- 
fects govern the reversal. The topology of the magnetic helic- 
ity domains was investigated under three different conditions: 
(i) emergence of the helix (zero- field cooling through T2), 
and reversal of the existing helix under (ii) equilibrium con- 
ditions (quasi-static electric-field-induced reversal) and (iii) 
non-equilibrium conditions (dynamic reversal). 

We now have to compare the switching here observed in a 
multiferroic with that of ferroelectric or magnetic domain re- 
versal by the adjunct electric or magnetic field, respectively, 
in non-multiferroics. This will highlight some of the principal 
differences between the magnetoelectric switching processes 
in this work and the well-documented switching of, e.g., ferro- 
electric domains by electric fields. There is no explicit theory 
on time-resolved magnetoelectric domain reversal so that we 
take a first qualitative approach here. 

Models for the reversal of ferroelectric domains involve nu- 
cleation at defects, fast growth along the field direction, and 
slow growth perpendicular to the electric fieldP^The "fast for- 
ward growth" occurs at the speed of sound and is related to the 
local charge at the domain walls perpendicular to the sponta- 
neous polarization. In Fig. [4] the forward growth is absent. 
Domain walls parallel and perpendicular to the poling field 
propagate at about 1 m/s and, hence, at a speed three orders 
of magnitude below that of sound. We suggest that this is 
related to the smallness of the spontaneous polarization in the 
magnetically-induced ferroelectric pLL^I and the resulting dom- 
inance of other mechanisms determining domain- wall propa- 
gation. One such mechanism is the coupling of the sponta- 
neous polarization to the spin system, and Fig. [2] revealed that 
this coupling is rigid. Therefore, the energies competing in the 
magnetoelectric domain reversal are the electric-field energy 
£ and the magnetic anisotropy energy K. Here an orders- 
of-magnitudes estimate reveals KjS « 100 which corrobo- 
rates the inertness of the magnetic system for the electric field 
acting on it: The electric-field energy is orders of magnitude 
weaker than the energy barrier separating the states with op- 
posite magnetic helicity. 

The resulting magnetoelectric nature of the switching pro- 



cess is reflected by the manifestation of the magnetic easy 
direction in the dynamic domain structure in Figs. [2jd)-(f) 
which is much more pronounced than in Fig.[T] Measurement 
of the electric-field dependence of the domain reversal beyond 
the saturation field, including a verification of the validity of 
Merz's lawP could yield further insight but this is restricted 
by the low breakthrough voltage of the helium gas cooling the 
sample. 

V. SUMMARY AND CONCLUSION 

In summary, time-resolved SHG experiments on multifer- 
roic MnW04 crystals revealed that the electric-field-induced 
reversal of the magnetization in multiferroics with magnet- 
ically induced improper ferroelectricity is inherently not ul- 
trafast. On the one hand, rapid reversal is inhibited by the 
inertness of the magnetic system in the electric field which 
is expressed by the smallness of the electric-field energy § 
in comparison to the magnetic anisotropy energy. On the 
other hand, the ex pecta tion that fast electronic contributions 
to the polarization 24 25 might transiently decouple from the 
slowly reversing magnetization was not confirmed: Accord- 
ing to our experiments the ferroelectric polarization remains 
rigidly clamped to the magnetization throughout the switch- 
ing process. 

With respect to applications, for magnetoelectrically con- 
trolled giant-magnetoresistance sensors 28 the electric-field in- 
duced reversal of the magnetic order parameter does not have 
to be faster than observed here. For these applications Figs. [I] 
to [4] show that memory effects result in a highly reproducible 
magnetoelectric switching process. No signs of fatigue were 
observed in our experiments. High durability may therefore 
be a major advantage of magnetoelectric devices. In addition, 
understanding and control of the pinning mechanisms caus- 
ing the memory effect will be the key for tailoring the mag- 
netically induced ferroelectric s according to technological re- 
quirements. 

Finally, there is a clear need for a theory modelling the 
non-equilibrium dynamics of the voltage-induced reversal of 
magnetic domains in multiferroics like MnW04. Here our 
orders-of-magnitude estimate of the energy scales involved in 
the switching process provides a clue about the approach such 
a theory might take. 
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